Kinetic magnetoelectric effect in a 2D semiconductor strip due to 
boundary-confinement induced spin-orbit coupling 



O 
O 



oo 

(N 



Yongjin Jiang 1 and Liangbin Hu 2 

department of physics, Zhejiang Normal University, Jinhua, Zhejiang 321004, P. R- China 
Department of physics and Laboratory of photonic information technology, 
South China Normal University, Guangdong 510631, P. R. China 

In a thin strip of a two-dimensional semiconductor electronic system, spin-orbit coupling may be 
induced near both edges of the strip due to the substantial spatial variation of the confining potential 
in the boundary regions. In this paper we show that, in the presence of boundary-confinement 
induced spin-orbit coupling, a longitudinal charge current circulating through a 2D semiconductor 
strip may cause strong non-equilibrium spin accumulation near both edges of the strip. The spins 
will be polarized along the normal of the 2DEG plane but in opposite directions at both edges of the 
strip. This phenomenon is essentially a kinetic magnetoelectric effect from the theoretical points of 
view, but it manifests in a very similar form as was conceived in a spin Hall effect. 
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I. INTRODUCTION 



There has recently been much interest in a fascinating 
topic in the research community, namely spin Hall effect 
( SHE ). SHE is such a phenomenon that a transverse 
spin current is generated when a longitudinal charge cur- 
rent circulates through a sample, and if the sample has a 
thin strip geometry, the transverse spin current will cause 
non-equilibrium spin accumulation at both edges of the 
sampl e 1 ! 2 ! 3 . Such a phenomenon would be much useful 
in the context of semiconductor spin-based electronics ( 
spintronics ) because it might provide an effective way 
for generating spin currents and/or non-equilibrium spin 
density in a nonmagnetic semiconductors without use of 
ferromagnetic metals or ferromagnetic semiconductors, a 
principal challenge in semiconductor spintronics^. From 
the theoretical points of view, such a phenomenon can 
arise from either intrinsic spin-orbit ( SO ) coupling ( 
i.e., spin-orbit splitting of the band structure or ex- 
trinsic SO coupling ( i.e., spin-orbit dependent impurity 
scatterings in a semiconducting material, and cor- 
respondingly, the phenomena due to intrinsic SO cou- 
pling was termed intrinsic spin Hall effect and the phe- 
nomena due to spin-orbit dependent impurity scatterings 
termed extrinsic spin Hall effect. From the standpoint of 
spintronic applications, intrinsic SHE are more attractive 
since it is an intrinsic property of a semiconducting ma- 
terial and does not rely on spin-orbit dependent impurity 
scatterings 2,3 . Recently two experiments were reported 
on the observation of SHE. One is on n-doped bulk GaAs 
by Kato et ali and the other is on two-dimensional p- 
doped GaAs by Wunderlich et &\&. The phenomenon 
observed in two-dimensional p-doped GaAs by Wunder- 
lich et al. was believed to be an intrinsic spin Hall effect 
since the edge spin accumulation measured in a thin strip 
of such a system is insensitive to impurity scatterings 
in the weak impurity scattering regime 8 , an important 
feature of intrinsic spin Hall effect 2 * 3 -. In contrast, the 
phenomenon observed in n-doped bulk GaAs by Kato et 



al. was believed to have an extrinsic origin ( i.e., due 
to spin-orbit dependent impurity scatterings ) since the 
edge spin accumulation measured in a thin strip of such 
a system is several order of magnitude smaller compared 
with the theoretical predictions of Refs.[2-3]. Although 
substantial progresses have been achieved on the detec- 
tion of spin Hall effect, it should be noted that from the 
theoretical points of view there are still intensive debates 
on whether intrinsic spin Hall effect does can survive in 
a spin-orbit coupled system. For example, several re- 
cent theoretical works have argued that, except for the 
ballistic transport limit, intrinsic spin Hall effect cannot 
survive in a diffusive two-dimensional electron gas with 
Rashba spin-orbit coupling even in the weak impurity 
scattering h m Hj2ii2iiiii2ii&±!ii£. As to the physical un- 
derstanding of the recent experimental results^, some 
significant controversies also exist which needs further 
clarifications. For details please refer to Refs. [16-25]. 

In this paper, we investigate theoretically another 
kind of electric-field driven edge spin-accumulation which 
might occur in a thin strip of a two-dimensional electronic 
( electron or hole ) system. From the theoretical points 
of view, this phenomenon is essentially a kinetic magne- 
toelectric effect^ due to boundary-confinement induced 
spin-orbit coupling ( which will be called edge SO cou- 
pling below ). But very interestingly, it would manifest 
in a very similar manner as was conceived in an spin Hall 
effect. For example, the electric-field driven edge spin ac- 
cumulation in a thin 2DEG strip due to this phenomenon 
will be polarized perpendicular to the 2DEG plane but 
along opposite directions at both edges of the strip; and 
in the weak impurity scattering regime ( below a certain 
disorder strength ) the edge spin accumulation will not 
decrease as the disorder strength increases, thus it can 
survive even in the diffusive transport regime. These fea- 
tures are very similar to the recently discovered spin Hall 
effect, but the mechanisms involved in this phenomenon 
are very different from that of the usual spin Hall effect 
from the theoretical points of view. The results obtained 
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in the present paper might provide some new implica- 
tions to the proper physical understanding of the recent 
experimental results^. 

The paper is organized as following: in section II 
an edge SO coupling model describing the boundary- 
confinement induced SO coupling in a thin 2DEG strip 
will be introduced and some details of the theoretical for- 
malism used in the paper will be briefly explained, and in 
section III some numerical results will be presented and 
discussed. 



II. MODEL AND THEORETICAL 
FORMULATION 

The system considered in the present paper consists 
of a thin 2DEG strip connected to two ideal leads, as 
was shown in Fig. 1(a). According to the theory of the 
relativistic quantum mechanics, if the movement of an 
electron is confined by a spatially non-uniform potential, 
the spin and orbital degree of freedom of the electron will 
be coupled together—. The longitudinal direction of the 
strip will be defined as the x direction and the normal of 
the 2DEG plane defined as the z direction. For the sake 
of simplification, we assume that the confining potential 
is spatially uniform along the longitudinal direction of 
the strip, i.e., the spatial variation of the confining po- 
tential occurs only in the transverse direction of the strip 
( defined as the y direction ). Under this assumption, the 
SO coupling due to the transverse spatial variation of the 
confining potential will take the following forr 
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where a = (a x , a yi & z ) are the usual Pauli matrices, k is 
the wave vector of electrons, V denotes the usual gra- 
dient operator and V(y) is the transverse confining po- 
tential, which depends only on the y coordinates. From 
the effective Hamiltonian (1) one can see that, for an 
electron moving along the +x ( or —x ) direction, the 
SO coupling will tend to force the electron spin to align 
along the +z ( or — z ) direction, depending on the sign 
of dV(y)/dy. In the absence of external electric field, 
the numbers of electrons moving along the +x and —x 
directions are equal and no net spin density can be re- 
sulted. However, if an external electric field is applied 
in the longitudinal direction of the strip, the numbers 
of electrons moving to —x direction will be larger than 
the number of electrons moving to the +x direction, and 
hence a net spin density polarized along the +z or ( — z ) 
direction might be resulted from the SO coupling. Such 
an effect was dubbed the kinetic magnetoelectric effect 
in the literature^. For a 2DEG strip, if the strip width 
is much larger than the lattice constant, one can further 
assume that the confining potential V(y) varies substan- 
tially only in a narrow boundary region as was shown in 



Fig. 1 (b) . ( We assume that the two boundaries of the srip 
are located at y = ±5- ). Then the SO coupling exists 
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FIG. 1: (a) A 2D semiconductor strip connected to two ideal 
leads; (b) Profile of the transverse spatial variation of the con- 
fining potential in the 2DEG strip; (c) Profile of the trans- 
verse spatial variation of the spin-orbit coupling strength in 
the 2DEG strip. 



only near both edges of the strip, and due to the symme- 
try of the confining potential, the SO coupling coefficient 
( proportional to dV(y)/dy ) has opposite signs at both 
edges of the strip, as was shown in Fig. 1(c). Thus, the 
edge SO coupling model described by the effective Hamil- 
tonian (1) implies that the electric-field driven spin ac- 
cumulation will be polarized along opposite directions ( 
normal to the 2DEG plane ) at both edges of the strip, 
similar to a spin Hall effect. This will be confirmed by 
detailed numerical calculations below. In order to make 
a comparison with the widely studied intrinsic spin Hall 
effect in a Rashba two-dimensional electron gas, we can 
also include a Rashba SO coupling in our edge SO cou- 
pling model. Then the total Hamiltonian of the 2DEG 
strip will be 
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where the second term is the Rashba SO coupling ( which 
arises from the inversion asymmetry of the trapping well 
along the normal of the 2DEG plane ) and a the Rashba 
SO coupling constant. ( z is a unit vector along the nor- 
mal of the 2DEG plane ). A direct use of the Hamiltonian 
(2) is not convenient if one wants to take into account the 
effects of impurity scatterings properly, which are very 
important in the diffusive transport regime. So in our 
calculations we will transform the Hamiltonian (2) into a 
discrete form. The discrete version of the effective Hamil- 
tonian for the total system ( including both the leads and 
the 2DEG strip ) will read: 
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Here t — 2 m*a 2 IS ^ ne hopping matrix clement be- 
tween two nearest-neighbour sites and a the lattice con- 
stant in the 2DEG strip. ^ Ri — {C Ri ^, C' Rii i) is the 
annihilation operators of electrons at the lattice site Ri 
in the strip. C Pja is the annihilation operator of elec- 
trons with spin a at the lattice site pj in the lead p ( 
p = 1,2 ). (p n ,R n ) stands for a nearest-neighbor pair 
of lattice sites across the interfaces between the 2DEG 
strip and the leads. w Ri is the on-site energy in the 
2DEG strip. In a clean system without disorder, one 
usually sets w Ri to be zero. t R = a /2a is the Rashba 
SO coupling coefficient in the 2DEG strip, which is as- 
sumed to be site-independent. The last term in Eq. @ 
stands for the boundary-confinement induced SO cou- 
pling and tft. = 4 ^ 2( ,; Vy is the coupling coefficient, 
which is site-dependent. The site-dependence of tjj will 
be determined by the actual form of the transverse con- 
fining potential in the 2DEG strip. In our calculations we 
assume a parabolic confining potential near both edges 
of the strip, as was shown in Fig. 1(b). In such cases, the 
site-dependence of the edge SO coupling coefficient can 
be expressed as 

t*. = ±t B max(ATg + \y\ - N y /2, 0), (4) 

where N y = W/a is the width of the strip (in units of 
lattice constant) , N B is the width of the narrow bound- 
ary regions in which the edge SO coupling exists, and 
t B is the minimum value of the site-dependent edge SO 
coupling coefficient. The sign ± is different for the two 
edges, as shown in Fig. 1(c). 

Our calculations will be based on the Landauer- 
Buttiker formula. To this end, we first consider the trans- 
mission and reflection of an electron incident from a lead. 
The real space wave function of an incident electron with 
spin a will be denoted as e~ tkmXp Xma(yp)> where XmaiVp) 
denotes the m'th transverse mode with spin index a in 
the lead p and fc^ the longitudinal wave vector. We adopt 
the local coordinate scheme for all leads. In the local co- 
ordinate scheme, the longitudinal coordinate x q in the 
lead q will take the integer numbers 1, 2,...,oo away from 
the 2DEG interface and the transverse coordinate y q take 
the value of — N q /2, N q /2. The longitudinal wave vec- 
tor kf n satisfys the relation — 2t cos(fc^) = E, where 
is the eigen-energy of the m'th transverse mode in the 
lead p and E the energy of the incident electron. Includ- 
ing both the incident and reflected waves, the total wave 
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function in the lead q has the the following general form: 
V<r{x q ,y q ) = 8 pq 5 aa ,e-^x p m M 

+ E<C> iKx *X Q nM (5) 

where (p P q ™a' stands for the scattering amplitude from the 
(ma) mode in the lead p to the (no~') mode in the lead 
q. To obtain the scattering amplitudes <p v ™°i from the 
Schrodinger equation ( which has now a lattice form 
and hence there is a separate equation for each lattice 
site and spin index ), we must solve the wave function 
ip^^iRi) in the 2DEG strip simultaneously. As deal- 
ing with usual scattering problems in quantum mechan- 
ics, we use boundary conditions to determine the scat- 
tering amplitudes (pg™^- In the lattice formalism, the 
wave functions in the entire 2DEG strip and in the first 
row of the leads ( i.e., x q = 1 ) will be involved in the 
boundary conditions, which will be determined from the 
Schrodinger equation. Because Eq.JSJ is a linear com- 
bination of all out-going modes with the same energy 
E, the Schrodinger equation is satisfied automatically in 
the lead q, except for the lattice sites in the first row ( 
i.e., x q = 1 ) of the lead which are connected directly 
to the 2DEG strip. The wave function in the first row 
of a lead ( which are determined by the scattering am- 
plitudes <fi q ™a' ) mus t be solved simultaneously with the 
wave function in the 2DEG strip. To simplify the nota- 
tions, the wave function in the 2DEG strip will be defined 
as a column vector ip whose dimension is 2N ( N is the 
total number of lattice sites in the 2DEG strip ). The 
scattering amplitudes 4^a' wm ^ e arranged as a column 
vector cf> whose dimension is 2M ( M is the total number 
of lattice sites in the first row of the leads ). From the 
Schrodinger equation for the 2DEG strip and for the first 
row of a lead, one can show that the boundary conditions 
at the interface between a lead and the 2DEG strip can 
be expressed as the following forms: 

Aip = h + B(p, 

C<p = d + D?/>, (6) 

where A and C are two square matrices with a dimen- 
sion of 2N x 2N and 2M x 2M, respectively; B and D 
are two rectangular matrices describing the hopping in- 
teraction between the leads and the 2DEG strip, whose 
matrix elements will depend on the actual form of the 
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geometry of the system. The vectors b and d describe 
the contributions from the incident waves. The details of 
these matrices and vectors will be given elsewhere^. 

After obtaining all scattering amplitudes 4> P ^^' from 
the boundary conditions, we can calculate the charge cur- 
rent in each lead through the Landauer-Buttiker formula, 
Ip = { e2 /h) J2 q X^ 1>CT2 {Tpol Vq ~ Tg;V p ), where V q = 
fi q /(—e) is the voltage applied in the lead q and /i q is the 
chemical potential in the lead q, T^, are the transmission 
probabilities defined by T^, = Em.nlO'V/V 
and v prn = 2£sin(fc£j is the velocity for the m'th mode 
in the lead p. The spin current in each lead can be calcu- 



lated similarly, /; = -(e/4ir) £ g ^ [{Tg' -Tg>)V 9 - 
{TP£ 2 - T£f 2 )Vp]. With the wave function Tp p ™ a {Rj) in 
the 2DEG strip at hand, the non-equilibrium spin den- 
sity in the 2DEG strip can also be calculated readily by 
taking proper ensemble average, and the results can be 
expressed as 

{s a {Ri)) = ^- ^ P /y Pm ^2'(P p a m ' T * (Ri)^^" 10 '(R%)> 

pma a, {3 

(7) 

where (S a (Ri)) denotes the spin density at the lattice 
site R, in the 2DEG strip. 



III. RESULTS AND DISCUSSIONS 

In our calculations we will take the typical values of 
the electron effective mass m = 0.04m e and the lattice 
constant a = 3nrr t — . The chemical potentials in the 
leads will be set by fixing the longitudinal charge current 
density to the experimental value ( w 100/iv4/1.5/i77i ) as 
reported in Ref. [8] . The Fermi energy of the 2DEG strip 
will be set to Ef = —3.8t throughout the calculations. 

In Fig. 2(a) we show the typical pattern of the spatial 
distribution of the electric-field induced nonequilibrium 
spin density < S z > in the 2DEG strip obtained based 
on the edge SO coupling model. ( Other components 
of the spin density is zero in the case of the edge SO 
coupling model, which were not shown in the figure ). 
We have used 180 x 60 lattice sites in total ( includ- 
ing both the leads and the 2DEG strip ) in the calcula- 
tions and the 2DEG strip contains 60 x 60 lattice sites. 
From Fig. 2(a) one can see that, the spatial distribution 
of the electric-field induced nonequilibrium spin density 
in the 2DEG strip obtained based on the edge SO cou- 
pling model manifests in a very similar form as was con- 
ceived in a spin Hall effecii£&, i.e., the spins are polarized 
perpendicular to the 2DEG plane but along opposite di- 
rections on both edge of the strip and the spin density 
has two opposite extrema near both edges. By fixing the 
longitudinal charge current density to the experimental 
value ( k 100 fiA/1, 5 [im ) as reported in Ref. [8] and set- 
ting the edge SO coupling coefficient to t B ~ 0.03t, we 
found that the spin polarization in the 2DEG strip ob- 
tained based on the edge SO coupling model ( « 1% ) has 
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FIG. 2: The two-dimensional real space configuration of the 
electric-field induced non-equilibrium spin density (S z (x,y)} 
(in units of ^) in the 2DEG strip obtained based on (a) the 
edge edge SO coupling model and (b) the Rashba SO coupling 
model, resepctively. In Fig. 1(a) the Rashba SO coupling co- 
efficient tn is set to zero and the edge SO coupling coefficient 
t B is set to be t B = 0.03t. The width of the boundary regions 
in which the edge SO coupling exists is set to be N B = 5. 
In Fig. 1(b) the Rashba SO coupling coefficient tn is set to be 
tn — 0.15t and the edge SO coupling coefficient t B is set to be 
zero. In both cases, the 2DEG strip contains 60 x 60 lattice 
sites. 



roughly the same order of magnitude as the correspond- 
ing experimental values reported in Ref. [8]. To make a 
comparison with the widely studied Rashba SO coupling 
models, in Fig. 2(b) we have also plotted the spatial dis- 
tribution of the electric-field induced nonequilibrium spin 
density < S x > in the 2DEG strip obtained based on the 
usual Rashba SO coupling model. From Fig. 2(a) and 
2(b) one can see that the typical patterns of the spatial 
distribution of the electric-field induced nonequilibrium 
spin density < S z > are very similar in both cases. But 
it should be pointed out that, in the case of our edge SO 
coupling model, only the z component of the spin den- 
sity is non zero. In contrast, for the case of the usual 
Rashba model, all three components of the spin density 
are nonzero. ( In Fig. 2(b) we have plotted only the spa- 
tial distribution of the z component of the spin density for 
comparison ) . Another slight difference that can be seen 
from Fig. 2(a) and 2(b) is that, the transverse spatial dis- 
tribution of the spin density has a very regularly striped 
pattern in the case of the edge SO coupling model, but 
for the case of the usual Rashba model, the spin-density 
pattern is not much regularly striped. This slight differ- 



5 



ence arises from the fact that, in the case of our edge 
SO coupling model, the SO coupling exists only in a nar- 
row boundary region near both edge of the 2DEG strip 
( see the illustration shown in Fig.l(b-c) ) and was as- 
sumed to be uniform along the longitudinal direction of 
the strip. These assumptions lead to a regularly striped 
spin-density pattern as shown in Fig. 2(a). For the case 
of the usual Rashba model, the SO coupling exists in en- 
tire the strip ( i.e., the SO coupling coefficient is nonzero 
everywhere in the strip ), thus the spin-density pattern 
in the strip is not much regularly striped. 



L=20 



(b) 



L=40 



t B /t=0.03 V 



(d) 



L=120 



(71=0.03 



nt=o.o3 



-? MM 




10 20 30 



FIG. 3: The transverse spatial distribution of the longitudi- 
nally averaged spin density {S z (y)} in the case of the edge SO 
coupling model. The lattice sizes of the 2DEG strip are cho- 
sen as: (a) 20 x 20, (b) 40 x 40, (c) 80 x 80, and (d) 120 x 120. 
The edge SO couling coefheent is set to be t B = 0.03. Fig. 3(e) 
shows the dependences of the transverse spatial distribution 
of {S z {y)) on the edge SO coupling strength, where the lattice 
size is fixed to be 60 x 60. In all these calculations, the width 
of the boundary region in which the edge SO coupling exists 
is set to be N B = 5. 

Next, we study the dependence of the edge spin accu- 
mulation on the strip width and the SO coupling strength 
in the case of the edge SO coupling model. Because 
the spatial distribution of the spin density has a regu- 
larly striped pattern along the longitudinal direction, we 
can use an averaged value of < S z {x,y) >, defined by 

< S z (y) >= -£ L S z (x,y)dx, as a measure of the spin 
accumulation. In Fig.3(a-d) we plot the profiles of the 
transverse spatial distribution of the spin accumulation 
obtained in several different cases with different lattice 



sizes. From these figures one can see that, < S z (y) > 
oscillates inside the strip and has opposite signs on both 
edges of the strip, and the magnitude of the spin density 
will reach a maximum value ( denoted as S%j ax below 
) near both edges of the strip. As the strip width in- 
creases, the transverse spatial distribution of the spin ac- 
cumulation become sharper and sharper near both edges 
of the strip and the oscillations of the amplitude of the 
spin density inside the strip tend to be smeared, i.e., 
the spin accumulation will be localized near both edges 
of the strip if the strip width is much larger than the 
lattice constant. From these figures one can also note 
that, the order of the magnitude of the spin density near 
both edges of the strip remains unchanged as the strip 
width increases, suggesting that the electric-field induced 
edge spin accumulation due to boundary-confinement in- 
duced edge SO coupling can survive even in the diffusive 
transport regime, similar to the phenomenon reported in 
Ref. [8] . ( For a 2D semiconductor strip with Rashba spin- 
orbit coupling, it was generally believed that the intrinsic 
spin Hall effect can not survive in the diffusive transport 
regime^iiSiiiiiSiiiiii^ ). The dependence of the trans- 
verse spatial distribution of the spin density on the edge 
SO coupling strength was shown in Fig. 3(e), from which 
one can see that the profiles of the transverse spatial dis- 
tribution of the spin density do not change substantially 
as the edge SO coupling strength varies. 

Finally, we discuss the effects of random impurity scat- 
terings on the electric-field induced nonequilibrium spin 
density in the case of the edge SO coupling model. To 
include properly the effects of random impurity scatter- 
ings, we assume that the on-site energy uj # 4 in the 2DEG 
strip are randomly but uniformly distributed in an en- 
ergy region [— Wd,Wd], where Wr> is the amplitude of 
the on-site energy fluctuations, which characterizes the 
disorder strengthi2*i£. We will calculate the spin den- 
sity for a number of random impurity configurations and 
then do impurity average. In Fig.4(a-b) we show the rela- 
tion between the edge spin accumulation and the disorder 
strength in both cases of (a) the edge SO coupling model 
and (b) the usual Rashba model, respectively. We have 
done impurity average over 10000 random impurity con- 
figurations for each case. From Fig. 4(a) one can see that, 
for the case of the edge SO coupling model, the edge spin 
accumulation does not decrease ( or may even increase ) 
as the disorder strength increases in the weak impurity 
scattering regime ( i.e., below a certain disorder strength 
) and for a fixed longitudinal change current density ( 
~ 100/1^4/ 1.5/xm ), similar to the intrinsic spin Hall ef- 
fect observed by Wunderlich et &\&. In contrast, for the 
case of the Rashba model, the edge spin accumulation de- 
creases monotonously as the disorder strength increases 
even in the weak impurity scattering limit, which can 
be seen clearly from Fig. 4(b). The different behaviors in 
these two models can be understood qualitatively as fol- 
lowing. For the edge SO coupling model, the SO coupling 
exists only in a narrow boundary region near both edges 
of the strip ( see the illustration shown in Fig.l(b-c) ). 
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FIG. 4: The dependences of the edge spin accumulation 
i^Max) on the disorder strength Wd in both cases of (a) the 
edge SO coupling model and (b) the Rashba SO model. In 
the weak impurity scattering regime ( below a certain disor- 
der strength ), the edge spin accumulation do not decrease 
as the disorder strength increases in the case of the edge SO 
coupling model. 



the strip will affect substantially the spin density. In con- 
trast, for the case of the Rashba model, the electric-field 
induced nonequilibrium spin density is due to the flow of 
a transverse spin Hall current, which will be damped 
significantly by all random impurity scattering events oc- 
curred in entire the strip and hence the edge spin accu- 
mulation will be decreased substantially with increasing 
disorder strength even in the weak impurity scattering 
limit. Of course, because localization effects will become 
important in the presence of strong impurity scatterings, 
the electric-field induced nonequilibrium spin density in 
the case of the edge SO coupling model will also be de- 
creased substantially in the presence of strong impurity 
scatterings, which can be seen clearly from Fig. 4(a). 

In summary, we have presented a microscopic model 
calculation for the kinetic magnetoelectric effect in a 
thin strip of a two-dimensional electronic system due to 
boundary-confinement induced edge SO coupling. We 
have shown that this effect can manifest in a very similar 
form as was conceived in a spin Hall effecti^, and some 
important features of this effect are similar to the in- 
trinsic spin Hall effect observed recently in thin strips of 
two-dimensional p-doped semiconductors^. The results 
obtained in the present paper may provide some new im- 
plications to the proper physical understanding of the 
recent experimental results. 
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